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We present a robust scheme by which fractional quantum Hall states of bosons can be achieved 
for ultracold atomic gases. We describe a new form of optical flux lattice, suitable for commonly 
used atomic species with groundstate angular momentum J g — 1, for which the lowest energy band 
is topological and nearly dispersionless. Through exact diagonalization studies, we show that, even 
for moderate interactions, the many-body groundstates consist of bosonic fractional quantum Hall 
states, including the Laughlin state and the Moore-Read (Pfaffian) state. Under realistic conditions, 
these phases are shown to have energy gaps that are larger than temperature scales achievable in 
ultracold gases. 

PACS numbers: 



There is intense interest in finding new settings in 
which topological phases of matter analogous to frac- 
tional quantum Hall (FQH) states appear. Ultracold 
atomic gases are ideal systems with which to achieve this 
goal: they allow studies of strong correlation phenom- 
ena for both fermions and bosons, and FQH physics can 
be approached for homogeneous fluids [T] as well as for 
atoms confined in optical lattices [2J. 

While existing theories of FQH-like phases in lattices 
have focussed on tight-binding models [HH9], one of the 
most promising routes to topological flat bands for ultra- 
cold atoms is through optical flux lattices (OFLs) [TDHT2] . 
An OFL uses a set of laser beams to produce a spatially 
periodic atom-laser coupling that induces resonant tran- 
sitions between two (or more) internal atomic states. The 
resulting energy bands, in particular the lowest one, have 
non-zero Chern numbers, and can be made narrow in en- 
ergy |12j . This opens the path to experimental studies 
of novel strong correlation phenomena in topological flat 
bands, notably the FQH effect of bosons. 

We present in this paper the first characterization of 
the many-body ground state of bosons in an OFL. We 
start with the design of a novel type of OFL, which 
fully exploits the structure of the most commonly used 
(bosonic) atomic species. The scheme is robust since, by 
contrast to some other OFL proposals [10], it does not 
require phase locking between the various optical beams 
composing the lattice. For optimized parameters its low- 
est band has Chern number 1 and is nearly dispersionless, 
closely analogous to the lowest Landau level for charged 
particles moving in a uniform magnetic field. We use 
exact diagonalization to determine the many-body spec- 
trum of a bosonic gas in this OFL. We show that FQH 
ground states appear for relatively weak atom interaction 
at the same filling factors as for a continuum Landau level 
PQ. Our work provides a concrete experimental scheme 
by which FQH states of bosons can be realized with large 
energy scales. Furthermore, it provides the first example 
of a non-Abelian quantum Hall state (the v = 1 Moore- 



Read state [13]) in a lattice model at high particle density 
with only two-body interactions. 

We focus in this paper on the case of atoms whose 
internal ground level has angular momentum J g = 1. 
This is the case for several stable bosonic isotopes of alkali 
metal species, namely 7 Li, 23 Na, 39 K, 41 K, 87 Rb. Wc 
denote \X), \Y), \Z) a basis of this level, defined such that 
Jx\X) — (and similarly for Y and Z). Here the set of 
directions X, Y, Z represents an orthogonal trihedron of 
the physical space [see Figure [lj a)] and Jx stands for the 
component of the angular momentum operator along the 
X direction. Note that one can replace \X), \Y), \Z) by a 
triplet of internal states selected among a more complex 
level scheme, as proposed e.g. in |14j . Our scheme will 
apply as long as each pair of states can be coupled by a 
resonant two-photon Raman transition with a negligible 
spontaneous emission rate [15] . 

We assume that \X), \Y), \Z) are the eigenstates of 
the atomic Hamiltonian in the absence of atom-laser 
coupling. We suppose that these three states are non- 
degenerate and non-equally spaced, and their energies are 
such that E x < E Y < E z , with E z - E Y ^ E Y - E X - 
For alkali atoms this situation can be reached by illu- 
minating the atomic sample with microwaves close to 
the hyperfine resonance (see supplementary information) . 
We denote by z the (1, 1, 1) direction of the X, Y, Z tri- 
hedron, and assume that the centre-of-mass motion of 
the atoms along the direction z is frozen. Therefore we 
consider in the following only the atomic motion in the 
perpendicular xy plane [see Figure [lja)]. 

The atoms are irradiated with laser beams propagating 
in the xy plane along three directions making an angle 
of 2tt/3 with each other. The three wave vectors are 
ki = k/2 (v3 u x + u y ), k-2 = k/2 (— y3 u x + u y ) and 
k% = —ku y , where {u x , u y } is an orthogonal unit basis of 
the xy plane. Here k stands for the typical wave number 
of the laser beams [16 . We choose the frequency com- 
ponents in each laser beam so that an atom can undergo 
resonant Raman transitions between the three internal 
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FIG. 1: (a) Atoms with a ground state with angular momentum J g = 1 are irradiated by three laser running waves propagating 
in the xy plane, whose wave vectors ki, i = 1, 2, 3 make an angle of 27r/3 with each other, (b) Triplet of light frequencies oji en- 
suring that the three possible Raman transitions are resonantly driven, (c) Graphic representation of three internal+momentum 
eigenstates, which are resonantly coupled by the laser beams whose frequencies are shown in (b). (d) Infinite array of inter- 
nal+momentum eigenstates that are resonantly coupled when three triplets of frequencies uii (red), uj[ (green), u" (blue) are 
simultaneously applied (see also Table El. 



states, by absorbing a photon in one wave and emitting 
a photon in a stimulated manner in another wave. The 
momentum change in such a transition is ±Qj, where 
qi = ki — ki+i- Here we set k^ = k\ and take h = 1 
for simplicity. 

Suppose first that each beam i consists only of a 
monochromatic plane wave with frequency and phase 
<Pi, and that the Wj's are chosen such the three Ra- 
man conditions are fulfilled: uj± — ll>2 = Ey — Ex, 
U)2 — u>3 = Ez — Ey [and thus u>\ — u>3 = Ez — Ex, 
see Figure |TJb)] . Each family of momentum eigenstates 
T{p) = {\X,p- fei), \Y,p-k 2 ), \Z,p- fe 3 )} generates a 
manifold that is globally stable with respect to atom-laser 
coupling. The three states of F{p) form an equilateral 
triangle in momentum space [Figure [ljc)]. The coupling 
between the atom and the laser field can be written (see 
supplementary material) 

V=-U ( \Y)(X\e l{qi - r+Vl - V2) + \Z)(Y\e i(q2 - r+V2 - V3) 

+ \X)(Z\e^ q3 - r+ ' P3 ^' P1 ^ + H.c, (1) 

where H.c. stands for Hermitian conjugate. The ampli- 
tude and sign of the coupling strength fl can be adjusted 
by tuning the intensity of the coupling lasers, and their 
detuning with respect to the atomic resonance. The fact 
that all three Raman transitions in Eq. ([!]) have the 
same amplitude is ensured by (i) taking the same in- 
tensity for each laser beam, (ii) choosing in-plane linear 
polarizations. A similar ring-coupling scheme has been 
used in [T7j to implement the Peierls substitution in a 
ID optical lattice. However in [T7] only two laser Raman 
transitions were used and the ring was closed using radio 
frequency transitions, which is not appropriate for our 
purpose. 

With only one triplet of laser frequencies as in 
Fig. [ljb), we do not produce the desired infinite peri- 
odic lattice for the atomic motion in momentum space 
[12] . However this goal can be reached by adding in- 
side the beams i two other triplets of frequency compo- 
nents Cc^ and ui" , i = 1,2,3. Here the roles are circu- 
larly exchanged with respect to the first triplet Wj: The 
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TABLE I: Phases of the Raman coupling matrix elements. 
Each line corresponds to a given momentum kick qi — ki — 
fei+i, and each column to a given pair of internal atomic 
states. This 3x3 array can be understood as a determi- 
nant: each of the six terms appearing in the calculation of 
this determinant corresponds to one of the six types of tri- 
angles in Fig. [TJd). The terms with positive (resp. negative) 
sign in the determinant calculation are for the upwards (resp. 
downwards) pointing triangles. 



u/'i (resp. uj") are such that uj' 2 — uj' 3 = Ey — Ex and 
u}' 3 — = Ez — Ey (resp. uj 3 — uj'{ = Ey — Ex and 
to" — uj 2 ' = Ez — Ey). In the following we suppose that 
the differences between the average frequencies u>, 
of the triplets are much larger than the splittings E a —Ep. 
Processes involving the absorption of a photon from a fre- 
quency triplet and the stimulated emission of a photon 
in another triplet thus play a negligible role. 

With the three frequency triplets acting simultane- 
ously on an atom, the family of states that are coupled 
to a given initial state can be represented by the infinite 
lattice in momentum space shown in Fig. |TJd) . Since 
there are 3 possible Raman transitions and 3 possible 
pairs of beams to induce a given transition, the atom- 
laser coupling V generalizing |l]) is now characterized 
by 9 matrix elements. These elements depend on the 9 
phases ifi,^, tp" and are summarized in Table 1 . From 
this Table, it is straightforward to write down explicitly 
the coupling V . For example the three terms appearing 
in ([!]) correspond to the diagonal terms of the array of 
Table H 

In order to characterize the possible non-trivial topol- 
ogy associated with the lattice in momentum space, we 
now evaluate the total phase gained by an atom when 
it undergoes a series of Raman transitions X — > Y — > 
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Z — > X and performs a closed loop in momentum space. 
This corresponds to traveling around the three sides of 
one of the triangles of Fig. [T|d). The resulting phase 
is different for upwards pointing triangles [such as the 
one of Fig. [ljc)] and downwards pointing ones [like the 
triangles labelled a,/3,7 in Fig. [ljd)]. For an upwards 
pointing triangle, the global phase is always zero. In- 
deed moving around the sides of such a triangle involves 
absorption and stimulated emission of photons whose fre- 
quencies belong to the same triplet, e.g. the uj^ triplet 
for the triangle of Fig. |ljc). Therefore each laser phase 
(fi (or ipt, ip'l) enters both with a + and a — sign in the 
total accumulated phase around such a triangle, leading 
to a null result. 

Downwards pointing triangles on the other hand cor- 
respond to a non-trivial phase. Consider for example the 
clockwise oriented path around the sides of the triangle 
labelled a in Fig.Qd): (i) The X -> Y transition is ac- 
companied by a change of atomic momentum q 2 , and it 
corresponds to a phase change ip' 2 — <p' 3 (see Table [I]); (ii) 
the Y — > Z transition is along q\, with the phase change 
(pi — cp 2 ; (hi) the Z — > X transition is along q 3 , with the 
phase change </?3 — tp±. As a result, the phase accumulated 
when traveling around the sides of triangle a is 

$ a = <p"-(pi + <p' 2 -<P2+(p3-V3- (2) 

We can similarly calculate the phases 3>0, 7 for the two 
other downwards pointing triangles. Although the sum 
^a+^p+^j is always zero, we can identify configurations 
such that each of these three phases takes a non-trivial 
value. For example the choice <pi — 27r/3, ^3 = —2tt/3, 
and all other phases equal to zero yields 

$ Q = §p = $ 7 = 2tt/3 mod 2vr. (3) 

From now on, we will stick to this choice, together with 
the assumption that fi > 0, which is obtained for an 
alkali-metal atom by tuning the lasers between the D\ 
and D 2 resonance lines. 

The OFL formed in this way has a reciprocal lattice 
spanned by the basis vectors G\ = 3<?i and G 2 =92- The 
real space lattice vectors are ai — ^J^ q (v3 u x +u y ) and 

a 2 = ^/§g U 2/) where q = \qi\ = \/3k. This geometry is 
equivalent to that of the three-state triangular flux lattice 
of Ref. 12. However, the pattern of phases in the recipro- 
cal space tight binding model differs: here we have fluxes 
of and 27r/3 in the upwards and downwards pointing 
triangles, as opposed to ir/3 for each TJ]. Nevertheless, 
the physical properties of the OFLs are very similar: in 
each unit cell of the real space lattice the lowest energy 
dressed state experiences = 1 flux quantum; the re- 
sulting bandstructure shows low energy bands that are 
analogous to Landau levels. In particular, the lowest 
energy band has Chern number of 1, and very narrow 
energy width, W, over a broad range of lattice depths 
fi. Here we focus on a lattice of depth fi = 3_Er [where 



En = q 2 /(2m) is the recoil energy for atomic mass m] 
close to which this bandwidth has a (local) minimum of 
W ~ 0.015-Er . In view of this very small bandwidth, the 
system is highly susceptible to the formation of strongly 
correlated phases even for relatively weak interactions. 

We have used exact diagonalization to study the 
groundstates of interacting bosons occupying the low- 
est energy band of the OFL for fi = 3-Er. (We ne- 
glect the population of higher bands, since the gap to 
the next band is very large, A ~ 46 W .) We consider the 
bosons to interact via spin-independent contact interac- 
tions, which is a good approximation for 87 Rb. We write 
the two-dimensional coupling constant as <?2D = 5r<7) 
where g is dimensionless. For atoms with 3D scatter- 
ing length a s restricted to 2D by a harmonic confinement 
of oscillator length ao, and neglecting (sub)band mix- 
ing, this is g = ^/8na s /ao [IS]. We study a finite sys- 
tem in a periodic geometry, with sides L\ = N\a\ and 
L 2 = N 2 a 2 , where N\^ 2 are integers. The total flux is 
then — N\N 2: so for N particles the Landau level fill- 
ing factor is v = N/N^. The interacting many-particle 
states can be classified by a conserved crystal momen- 
tum at Ncf, points in the Brillouin zone. We construct 
the Hamiltonian at each crystal momentum and use a 
standard Lanczos method to determine the low energy 
spectrum. 

For very weak interactions, g -C 1, the bosons form a 
condensate in the minima of the band dispersion. How- 
ever, our numerical results show that, even for moder- 
ate interaction strength g > 0.2, this (compressible) con- 
densed phase is replaced by strongly correlated (incom- 
pressible) FQH states at filling factors v = 1/2,2/3,3/4 
and 1. Here, we focus on the FQH states at v = 1/2 
and 1. (Results for v — 2/3,3/4 are described in the 
Supplementary Material.) 

Evidence for the appearance of incompressible phases 
is found by calculating the discontinuity in the chem- 
ical potential A/x for the groundstate: the difference 
between the chemical potential for adding a parti- 
cle and that for removing a particle. A non-zero 
and positive Afi indicates that the system is incom- 
pressible. To minimize finite-size effects we define [TU] 
A/i = N[E N+1 /(N + 1)+E N _ 1 /(N-1)-2E N /N} 1 
where En is the groundstate energy for N particles. 

In Fig. [2] we plot the dependence of A/i on interaction 
strength g at filling factors v = 1/2 and 1. For v = 
1/2 there is an onset of incompressibility for g > 0.2, 
and for v = 1 incompressibility appears for g > 0.4. 
In the thermodynamic limit, N — > 00, the transitions 
from compressible Afi = to incompressible A/i > 
should be sharp, and can even be discontinuous for first- 
order transitions, but are rounded in Fig. [2] by finite-size 
effects. The observed rises of A/i are indications of the 
approximate values of g at which there are transitions 
into the incompressible phases. 

To explore the nature of these incompressible phases it 
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FIG. 2: Incompressibility, as measured by the discontinuity 
in the chemical potential Afi defined in the text, as a function 
of interaction strength g at several filling factors: v = 1/2 (for 
N = 9 bosons in a system of size N$ = N1N2 = 6x3, circles); 
and v = 1 (JV = 12 in N# — 6x2, triangles). The inset shows 
A/x as a function of filling factor for a series of interaction 
strengths g. 



is instructive to study their (neutral) excitation spectra 
in the strong-interaction limit g — > 00. These spectra, 
Fig. |3j show all the expected properties of the bosonic 
Laughlin (y = 1/2) and Moore- Read [y = 1) states. On 
this periodic geometry, these topologically ordered in- 
compressible phases should show groundstate degenera- 
cies (of 2 and 3 respectively) in the thermodynamic limit, 
separated by an energy gap from the remaining excita- 
tions. As shown in Fig. [3j even for these finite systems 
these groundstate degeneracies appear clearly. Results 
on other system sizes and geometries (not shown) are 
consistent with these results, confirming that these near 
degeneracies are robust features, not imposed by symme- 
tries, that characterize these topological phases. 

The FQH states that we find for the OFL (at v = 
1/2,2/3,3/4, 1) are the same as those found for contact 
interacting bosons in the continuum lowest Landau level 
(LLL) pp. We have established the equivalence of the 
phases of these two models by studying the evolution of 
the many-body spectrum for a series of Bloch wavefunc- 
tions that interpolate between those of the lowest band 
of the OFL and those of the LLL. To do so, we consider 
a fictitious atom with N s — 12 internal states, and repre- 
sent the LLL by the iV s = 12 triangular OFL of Ref. H2 
the lowest band of which has properties that are indis- 
tinguishable from those of the LLL for suitable coupling 
f2' ~ IOE'r [5D]. We place 9 additional internal states 
at the midpoints of the bonds of Figjljd), coupled to 
each other and to the original states X, Y, Z by bonds of 
strength 0,' and with 7r/12 flux through each new triangu- 
lar plaquette [21]. Choosing (0,0') = (3(1 - A), 10\)E R 
and varying A leads to smooth interpolation of the lowest 
energy band and the many-body spectrum, from those of 
the present model (A = 0) to those of the LLL (A = 1). 
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FIG. 3: Low-energy spectra for the OFL with Q/Er = 3 
in the strong-interaction limit g — > 00 at filling factors (a) 
v = 1/2 (N = 9 bosons in = NiN 2 = 6x3) and (b) v = 1 
(N = 12 bosons in N1N2 = 6 x 2). The crystal momentum 
k = a 1 Gi/Ni+a2G2/N2 is labelled by the index, i = I+Q1 + 
JViQ2 for qi = 0, . . . , Ni — 1 and 02 = 0, -/V2 — 1. The quasi- 
degenerate groundstates have the expected multiplicities and 
crystal momenta for the Laughlin state (y = 1/2), and the 
Moore-Read state [y = 1). 



In all cases [y = 1/2, 2/3, 3/4, 1) the energy gap remains 
open, showing that the phases of these two models are 
the same [351 Indeed, there is very little change in 
the spectrum, showing that this OFL (with £1 = 3-Er) is 
a very close representation of the LLL. For example, for 
the LLL the v = 1/2 state has zero interaction energy, 
as the two-body correlation function vanishes exactly at 
zero range. Here, the v = 1/2 state in the OFL has 
E N /N ~ 7 x ICT^-Ek [see Fig. Ufa)] showing that the 
zero-range two-body correlation function nearly vanishes. 

To summarize we have proposed a robust atom-laser 
configuration that can lead to FQH states of bosons in 
a well- accessible range of parameters. The robustness of 
the setup is ensured by the absence of need for a sta- 
bilization of the relative phases of different beams |24) . 
The only phase difference to be controlled is within each 
single-mode laser beam (<Pi,<p'i,<Pi) and can be set by 
acoustic-optic modulators driven by programmable func- 
tion generators. The lowest energy band is insensitive 
to fluctuations in the laser amplitudes around Q = 3-Er, 
its bandwidth increasing by less than 10~ 3 Er within the 
range Q/En = 2 — 4, and its topology remaining un- 
changed. The minimal interaction strength g s» 0.2 for 
obtaining FQH states corresponds to a 2D confinement 
frequency of > 7 kHz for Rb, which is readily achieved 
in an optical lattice. A clear signal of the formation of 
strongly correlated phases would be the appearance of 
density plateaus (wedding cake structure) in in-situ im- 
ages of the gas, arising from incompressibility A/i > 0. 
This requires the temperature to be smaller than A/i, 
which for the Laughlin state we find from Fig. [2] to be 
« 0.02£t> for g = 0.4, that is 10 nK for 87 Rb. 
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Supplementary material 



THE BASIS SET \X), \Y), \Z) 

The components Jx and Jy of a spin 1 angular mo- 
mentum operator read in the standard eigenbasis {\mz = 
+l),\m z = 0),\m z = -l}}oiJ z : 



Jx = 



V2 



/o 1 N 
1 1 

V° 1 °y 



Jy = 



V2 



-i 

1 -i 
,0 i 



(SI) 



In the standard eigenbasis \mz = 0, ±1), we define the 
eigenvectors |a) of J a (a — X, Y, Z) with zero eigenvalue: 





\Z) 



(S2) 

These three states also form an orthonormal basis set 
(denoted hereafter the Cartesian basis) and are such that 

J x \Y)=i\Z), J Y \Z)=i\X), J z \X)=i\Y). (S3) 

For any real 3- vector u, the state ux\X) +uy\Y) +uz\Z) 
is the eigenstate of u ■ J with eigenvalue 0. 



MICROWAVE DRESSING OF A SPIN-1 
GROUND STATE 

We consider an atom from the alkali-metal family with 
a nuclear spin I = 3/2, so that its ground level is split in 
two hyperfine levels with angular momentum F — 1 and 
F = 2. The atom is irradiated with a linearly polarized 
microwave (mw), whose frequency is detuned by A mw 
from the frequency tJhf of the F = 1 -H> F = 2 hyperfine 
transition. The coupling between the mw and the atom 
(half-Rabi frequency) is defined by K mw = /iB-3 mw /2, 
where fiB is the Bohr magneton and B mw is the ampli- 
tude of the oscillating mw magnetic field. We restrict to 
the case where n mw <C A mw <C wm- We can then treat 
the atom-mw coupling using the rotating wave approxi- 
mation and evaluate the shifts of the Zeeman states of the 
F = 1 level using second-order perturbation theory. Tak- 
ing the quantization axis parallel to the polarization axis 
(Z) of the microwave, the shifts of the states mz — 0, ±1 
read [T] 



AE(m z ) 



1 



(S4) 



The action of the mw on the F = 1 hyperfine level can 
thus be described by the operator a + f3J z , where j3 = 



Suppose now that a second mw, at a different detun- 
ing A^ w and with a linearly polarization along X, also 
irradiates the atom. In the perturbative framework used 
above, the combined action of the two mws is described 
by the effective Hamiltonian -ff m w,cff = PJ Z + P'Jx> U P 
to an additive constant. The eigenstates of i? mWje ff are 
the states \X), \Y),\Z) with the energies /?, (3 + /3', /?', 
respectively. 

To estimate the required value for B mwi we start 
by noticing that when writing the atom-laser coupling 
[Eqn (1) of the main text], we only take into account 
the resonant elements. For example we assume that the 
atom can undergo a transition \X) — > \Y) via the absorp- 
tion of a photon k\ and the emission of a photon k2, and 
we neglect the transitions \X) — > \Y) occurring via the 
absorption of &2 and emission of k\. This is legitimate 
when the two-photon coupling fi is small compared to 
the energy difference between \X), \Y), \Z): Q <§C /3,/3'. 
We also note that the perturbative approach leading to 
(S4) is valid only when K mw <§; A mw . Requiring that 
the two sides of each strong inequality differ by at least 
a factor 5, we find ft < K mw /100. Consider as an illus- 
tration the case of Rubidium atoms for which the opti- 
mal two-photon coupling ft = 3q 2 /(2m) = 9k 2 /(2m) « 
2ir x 30 kHz. The microwave coupling K mw has to be 

> 3 MHz, corresponding to a microwave magnetic field of 

> 4 G. This is a large, but still realistic value, especially 
if one uses small resonant loops to increase the value of 
B mw at the location of the atomic sample. 



THE LIGHT-SHIFT OPERATOR FOR AN 
ALKALI-METAL ATOMIC SCHEME 

We consider an atom that is irradiated by a monochro- 
matic laser beam of frequency U). The light-shift operator 
V gives the restriction of the atom-laser coupling to the 
ground atomic level g at first order in laser intensity [2]. 
In the absence of electron and nuclear spins, the ground 
and excited states that are involved in the resonant tran- 
sition of an alkali- metal atom have an angular momentum 
J g = and J e — 1, respectively. The atom-laser coupling 
can be written 



U = K a \e,a){g\+R.c. (S5) 

a=X,Y,Z 

= Yl «m|e,m)(c?|+H.c. (S6) 



m=0,±l 



Here the 



= X } Y, Z, (resp. K m , m = 0, ±1) denote 



the atom-laser coupling strengths (half-Rabi frequencies) 
in the Cartesian (resp. standard) basis with 

1 



K± 



V2 



K a = K Z - 



(S7) 



mw 

/(4A mw ). 



In this "no-spin" case, the light-shift operator simply de- 
scribes the displacement of the ground state by the quan- 
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FIG. SI: Amplitude of the couplings for the two components 
Di (top) and D2 (bottom) of the resonance line of an alkali 
atom. Here the spin of the nucleus is supposed to be zero. The 
coupling amplitudes for the g- state take symmetric values 
with respect to those indicated for the g + state. 



tity |k| 2 /A, where A is the detuning of the laser with 
respect to the atomic transition, k = (kj, Ky, kz) and 



E 



=X,Y,Z 



E 

m=0,±l 



(S8) 



We now take the spin of the valence electron into ac- 
count, but still assume that the nucleus has spin zero. 
The ground level is then a two-dimensional manifold 
(angular momentum J g — 1/2). Because of the fine- 
structure Hamiltonian, the resonance line of the atom is 
split into two components Di and D2, corresponding to 
the transition from the ground state to the excited states 
with total angular momentum J e = 1/2 and J e = 3/2, 
respectively. 

Let us focus first on the light-shift operator associated 
with the Di transition, which is dominant if the detuning 
Ai = bj — uj\ of the laser excitation from the Di line (fre- 
quency oji) is much smaller than the detuning A2 from 
the D2 line. Using the well-known Clebsch-Gordan coef- 



ficients (see e.g. Figure SI) and taking Z as quantization 
axis, we get in the corresponding {|<?+), \g~)} basis: 



y = 1 ( |k | 2 + 2|k_| 2 -V2(k k*_ +k+k,* ) 
3Ai {-V2(k^k^ + k qK *_) \k \ 2 + 2\k+\ 2 



This coupling can be written in a compact form 



(S9) 



(S10) 



where 1 and S are the identity and spin operators for the 
spin 1/2 ground state, respectively, and where the (real) 
effective field Bi for the H\ transition is: 

B 1 = -^-KXK*. (Sll) 



We now take into account both the Di and D2 transi- 
tions. A straightforward generalization of the preceding 
calculation leads to 



V = A1 + B-S, 



with 



and 



A = 



B = 



2i 



1 

AT 



1 

AT 



1 

AT 



2 

AT 



K X K 



(S12) 



(S13) 



(S14) 



As a final step we take into account the nucleus spin /. 
As above we consider the case I — 3/2, which leads to a 
splitting of the ground level into two hyperfme states of 
angular momentum F = 1 and F — 2. Here we consider 
the F = 1 state, with the three Zeeman sublevels 



\F = l,m F = ±l) 
\F = l,m F = 0) 



= T- 



V3 



1 3. 1, 1 1, 
2 1 2 2' 2 1 2 2' 



71' 



2'2 y 



where the state |m e ; m n ) is labelled by the quantum num- 
bers for the projection along the Z axis of the electron 
spin (m e ) and nuclear spin (m„). We suppose that the 
detunings Ai and A2 of the laser with respect to the ex- 
cited states J e = 1/2 and J e = 3/2 are large compared to 
the hyperfme splittings of these levels. The calculation of 
the restriction of the coupling (S12 ) to the F = 1 ground 
level then gives: 



V = Ai + B' ■ F, 



where 



B' = - 



4 6 lAi 



K X K 



(S15) 



(S16) 



Using ( S3 ) the vector part of the coupling U v , 



U - 

A 1 can be written in the Cartesian basis for the F = 1 
ground state: 

Koc = i (B' Z \Y)(X\ + B' X \Z)(Y\ + B' Y \X) (Z\) + H.c. 

(S17) 



HAMILTONIAN FOR RESONANT RAMAN 
TRANSITIONS 



Consider now a scheme such as the one of Fig. S2 
where an external field lifts the degeneracy between 
the three states of the Cartesian basis, and where two 
monochromatic light waves with the same amplitude 
(characterized by k > 0) and with frequency, phase, wave 




FIG. S2: Laser scheme providing a resonant Raman coupling 
between two sublevels of the F — 1 ground state. 



vector and polarizations LJi,ipi,ki,€i, i — 1,2 induce a 
resonant Raman coupling between \X) and \Y). 

The transition from \X) to \Y) occurs resonantly via 
the absorption of a photon in wave 1 and the stimu- 
lated emission of a photon in wave 2. Keeping only this 
resonant process, the releva nt co ntribution to the cross- 
product n x n* entering in (S16) is k 2 (eie^ fcl ' r+¥ ' 1 ^) x 



\e.2e° v " z ■ • 1-*/ j f where r is the position of the atom. 

In this work we restrict to waves with linear (real) po- 
larizations in the xy plane, making an angle of 27r/3 with 
each other. In this case e\ x e 2 — (y/3/2)u z . The projec- 
tion of u z on each basis vector of the trihedron X, Y, Z 
is l/v3, and the value of B' z that is relevant for the res- 



onant coupling between \X) and \Y) (see (S17l) is thus 

(S18) 



.2 



B' 



IK 

12 



1 



1 

A~2 



i{qx-r+ipx-ip 2 ) 



with q 1 = ki — fc 2 . This leads to the resonant part of the 
vector light-shift operator: 



V^fO = -n e i(n-r+vi-V2) \Y)(X\ +H.C. 



where we have set 

n 



K 

12 



1 



1 



(S19) 



(S20) 



The coupling strength Q is positive when the laser is 
tuned between the Di and the D 2 lines (Ai > > A 2 ), 
and negative otherwise. 

Consider for example the case of Rubidium atoms and 
choose the detuning such that the scalar part of the 
atom laser coupling (S13) vanishes: A 2 = — 2Ai, cor- 
responding to the wavelength A = 790 nm (Ai = 795 nm, 
A 2 = 780 nm). The optimal two-photon coupling is 
tt = 3£r = 9fc 2 /(2m) fa 2n x 30kHz and the photon 
scattering rate is 7 pa 9 x Tn 2 / (2 Af) ps 20 s -1 , where 
the factor 9 accounts for the number of monochromatic 
beams shining on the atoms. The time scale for establish- 
ing a many-body state such as the Laughlin state can be 
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FIG. S3: Low-energy spectra for the OFL with Q/E-r — 3 
in the strong-interaction limit g — > oo for a system of size 
N,/, — N1N2 — 6 x 2 at filling factors for which the groundstate 
is incompressible. The crystal momentum A; = aiGi/6 + 
a2G2/2 is labelled by the index, i = 1 + ai + 6Q2 for ai = 
0, ...,5 and Q2 = 0,1. The quasi-degenerate groundstates 
have the expected multiplicities and crystal momenta for the 
Laughlin state (y = 1/2), composite fermion states {y = 2/3, 
3/4), and the Moore-Read state [y = 1). 



estimated as the inverse of the corresponding gap A/x -1 . 
Taking A/i = 0.02 Er as a typical value (see Fig. 2 of the 
main text), we obtain A/i^ 1 ps 1ms. The heating due to 
photon scattering should then play a minor role during 
this time. 



ADDITIONAL NUMERICAL RESULTS 

We present in this section some additional numerical 
results alluded to in the main text. 

In Fig. S3 we present the excitation spectra for a 
system of size Ni x 7V 2 = 6 x 2 at all filling factors 



[y = 1/2,2/3,3/4 and 1) for which we find evidence for 
incompressible FQH states. These show the expected fea- 
tures of the Laughlin state (v = 1 /2) , composite fermion 
states [y = 2/3,3/4) and Moore-Read state (y — 1) of 
bosons. On this periodic geometry, these topologically 
ordered incompressible phases should show groundstate 
degeneracies (of 2,3,4 and 3 respectively) in the ther- 
modynamic limit, separated by an energy gap from the 
remaining excitations. This is indeed the case, as may 



be seen clearly in Fig. S3 



In Fig. S4 we present the equivalent excitation spec- 
tra for a lowest band formed from states in the lowest 
Landau level (LLL). As described in the main text, we 
represent the LLL by the N s = 12 OFL lattice of Ref.H 
the lowest band of which has properties that are indis- 
tinguishable from those of the LLL for a lattice coupling 
of SI' = 10-Er. This very close equivalence arises from 
the very rapid convergence of the properties of the in- 
state OFL lattice of Ref . [3] to those of a charged particle 
in a uniform magnetic field, the spatial fluctuations of 
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FIG. S4: Low-energy spectra for contact-interacting bosons 
in the LLL for a system of the same size and geometry as 
Fig. [S3](iV = 12), at filling factors v = 1/2,2/3,3/4,1 for 
which the groundstate is incompressible. The crystal momen- 
tum is labelled as for Fig. |S3| 
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FIG. S5: Low-energy spectrum for the OFL with Q/E-r = 3 
in the strong-interaction limit g — > oo for a system of size 
N<f, = iViiV2 = 7 x 2 at filling factor v = 1. The crystal 
momentum fc = aiGi/7 + aiGi!^ is labelled by the index, 
i = 1 + Qi + 7o?2 for a,\ — 0, . . . , 6 and 0:2 = 0,1. The 
three quasi-degenerate groundstates (marked by circles) are at 
the expected crystal momenta for the Moore-Read state. No 
symmetry protects these quasi-degenerate levels. The same 
quasi-degeneracy appears for contact interacting bosons in the 

LLL pg. 



the energy and effective magnetic field experienced by 
the lowest energy dressed state of the OFL falling expo- 
nentially with increasing N s , and already negligible for 
N s = 12 [3]. Interpolation between the N s = 3 OFL of 
this paper and the LLL leads to a smooth evolution both 
of the single-particle levels of the lowest bands, and of 
the many body spectrum for bosons occupying this low- 
est band. This establishes that the FQH phases of these 
two systems are the same. Moreover, the spectra of the 
OFL (Fig. S3 ) differ only very slightly from those of the 
The main qualitative difference is that 



S4| 



LLL (Fig. 

some approximate degeneracies in Fig. [S3] become exact 
degeneracies connected to the many-body translational 
symmetry of the LLL [3]. Furthermore, the approximate 
3-fold groundstate degeneracy of Fig. [S3^d) becomes an 
exact degeneracy in the LLL Fig. |S4[d). This is due to 
a 7r/3 rotational symmetry of the LLL in this geometry 
of Nx x A^ 2 = 6 x 2 (for which the sides of the simulation 
cell have equal length \Li\ = I-&2I), which transforms the 
three groundstates of the Moore-Read state in the LLL 



into each other. We emphasize that, in general, this 
three-fold quasi-degeneracy is unrelated to any symme- 
try, but is an emergent quasi-degeneracy in the thermo- 
dynamic limit. This is evidenced by studies on other 
system sizes and geometries, just as in the LLL [5]. For 
example, Fig. [S5] shows the spectrum for the OFL for 
iV = 14 particles in a system of size Ni x N2 — 7 x 2 for 
which no symmetry relates the three quasi-degenerate 
groundstates. 
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